Exciton-phonon information flow in the energy transfer process of photosynthetic complexes 
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Non-Markovian and non-equilibrium phonon effects are believed to be key ingredients in the energy transfer 
in photosynthetic complexes, especially in complexes which exhibit a regime of intermediate exciton-phonon 
coupling. In this work, we utilize a recently-developed measure for non-Markovianity to elucidate the exciton- 
phonon dynamics in terms of the information flow between electronic and vibrational degrees of freedom. We 
study the measure in the hierarchical equation of motion approach which captures strong system-bath coupling 
effects and non-equilibrium molecular reorganization. We propose an additional trace-distance measure for the 
information flow that could be extended to other master equations. We find that for a model dimer system and 
the Fenna-Matthews-Olson complex that non-Markovianity is significant under physiological conditions. 



The initial step in photosynthesis involves highly efficient 
excitonic transport of the energy captured from photons to a 
reaction center (T). In most higher plants and other organisms 
such as green sulphur and purple bacteria this process occurs 
in light-harvesting complexes which consist of electronically 
coupled chlorophyll molecules embedded in a solvent and a 
protein environment [2 |. Several recent experiments show that 
excitonic coherence can persist for several hundreds of fem- 
toseconds even at physiological temperature (3). These exper- 
iments suggest the hypothesis that quantum coherence may be 
biologically relevant for photosynthesis. 

The challenging regime of intermediate coupling of the 
electronic to the vibrational degrees of freedom motivates the 
study of sophisticated master equations with non-Markovian 
effects. Several approaches have been taken, ranging from a 
polaron transformation [4|, quantum state diffusion 0, gen- 
eralized Bloch-Redfield J6), non-Markovian quantum jumps 
HE), quantum path integrals (QUAPI) [9|, to density matrix 
renormalization group (DMRG) iflOl . After photoexcitation, 
the nuclear coordinates of a molecule will relax to new equi- 
librium positions, a phenomenon that is borne out by the time- 
dependent Stokes shift, that is the frequency shift between ab- 
sorption and emission spectra. Redfield theory assumes that 
the phonon bath is always in thermal equilibrium and thus can- 
not capture this effect ifTTI . Ishizaki and Fleming employed 
the hierarchical equation of motion (HEOM) approach lfl2l . 
which interpolates between the usual weak and strong (sin- 
gular) exciton-phonon coupling limits and takes into account 
non-equilibrium molecular reorganization effects. 

In this work, we study the non-Markovianity of the exciton 
transfer process by means of numerical simulation. Quantum 
mechanical time evolution under decoherence leads to trans- 
fer of information encoded in the electronic excited state to 
the vibrational degrees of freedom. Non-Markovianity can 
be characterized as the return of this information to the elec- 
tronic degrees of freedom. To quantify the exciton-phonon 
information flow, we employ the state-of-the art master equa- 
tion approach by Ishizaki and Fleming and a newly-developed 
measure for non-Markovianity. The central question we will 
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answer is how much information is exchanged between exci- 
tonic and phononic degrees of freedom in that process and, 
specifically, how much information returns from the phonons 
to the exciton. 

Measure for non-Markovianity— A recent area of research 
is the study of measures to characterize non-Markovianity ifPJl 
Q3). We utilize the readily applicable measure developed by 
Breuer and co-workers lfl3l . It is based on the quantum state 
trace distance: 

D( Pl (t),P2(t)) = \Tr{\ Pl (t)-p 2 (t)\}, (1) 

where \A\ — V A. This measure quantifies the distinguish- 
ablility of two quantum states Pl and p 2 . For a completely 
positive trace-preserving map E the trace distance is contrac- 
live, i.e. D(E( Pl ),E(p 2 )) < D( Pl ,p 2 ) [151. Time in- 
tervals in which D increases indicate non-Markovian infor- 
mation flow IT31 . in our case from the vibrational degrees of 
freedom back to the electronic degrees of freedom. In terms of 
the slope of D, i.e. a ( Pl (t) ,p 2 (t)) = f t D ( Pl (t) , p 2 (t)) , 
this implies that a > for these time intervals. Formally, 
a measure for non-Markovianity is defined by the following 
ED: 

NM-ity = max / dt 1 (a (t)) a ( Pl (t) , p 2 (i)) . 

Pi(0),P 2 (0) Jo 

(2) 

Here, the indicator function X (x > 0) = 1 (T (x < 0) = 0) 
monitors only the increasing part of the trace distance time 
evolution i.e. the information backflow. The optimization of 
the initial states in Eq. |2]) obtains the maximally possible non- 
Markovianity of a particular quantum evolution. In the case of 
photosynthetic energy transfer, initial states are given by the 
particular physical situation and lead to smaller values for the 
non-Markovianity. We use both physical and optimized initial 
states in this work. 

Hierarchy equation of motion— The Hamiltonian 
describing a single exciton in a complex with N 
molecules is given by H e — J2m=i( e ™ + X)\m){m\ + 
E m <„ J-mn (\m)(n\ + \n)(m\). The site energies e m , and 
couplings J mn are obtained from detailed quantum chemistry 
studies and/or fitting of experimental data. The set of states 
|m) denotes the site basis. The reorganization energy A is the 
energy difference of the non-equilibrium phonon state after 
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Franck-Condon excitation and the equilibrium phonon state 
and is assumed to be the same for each site. The superopera- 
tor corresponding to H e is C e p — [H e ,p]. The Hamiltonian 
for the phonon environment is iJ p h = Y^i \Pi + 9?) /2j 
where Pi, and ^ are the frequency, dimensionless mo- 
mentum and position operators of the phonon mode i. The 
coupling of the electronic degrees of freedom to the phonons 
is given by the Hamiltonian H cx - ph = J2 m v m {J2i 9i1i) m < 
with V m = \m)(m\ and where we assume that site-energy 
fluctuations dominate, the bath degrees of freedom at each site 
are uncorrelated, and the coupling constants to the respective 
modes are given by the constants gi. The spectral density, 
which describes the coupling strength of exciton to particular 
phonon modes, is given by J (cj) = 2A-fui/ (oj 2 + r ) 2 ). 
Perhaps the most relevant parameter, the bath dissipation rate 
7, is related to the bath correlation time by t c = 1/7. The 
equation of motion for the reduced single-exciton density 
matrix p is obtained by tracing out the less relevant phonon 
degrees of freedom and utilizing Wick's theorem for the 
Gaussian fluctuations, leading to a non-perturbative hierarchy 
of system and auxilliary density operators (ADOs) fl2l : 
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The system density matrix is the first member of the hierar- 
chy pit) — a (t) . The other members of the hierarchy are 

arranged in tiers and indexed by n = (m , njv) > 0, 

where a single tier is given by all elements where J2 m n ™. 
is constant. The notation n m ±i = (/ii n m ± 1, tin) 
is introduced. The superoperators in the high temperature 
regime fayf} < 1 are <j) m a n — i[V m ,a n ] and 8 m a n = 
i (2X//3H 2 [V m ,a n ] - iXy/h{V m ,a n }) and act on system 
and ADOs, where {, } is the anticommutator. Initially, all the 
hierarchy members are zero, which corresponds to the elec- 
tronic ground state phonon equilibrium configuration. For nu- 
merical propagation, the hierarchy is truncated based on the 
criterion J2 m Um ^> wnere w e is a characteristic fre- 
quency of C e . We assume that this condition is satisfied when 
the J2m n m ls by a factor 5 larger than ^ . 

The representation of the hierarchy equation of mo- 
tion in Eq. ^ is not unique. It can be advan- 
tageous to expand the hierarchy in a set of normal- 
ized auxiliary systems lfT6l . Redefining the ADOs as 

^{t) = (IL 

superoperators as (f> m a nm+1 = y/Jj 
and n m 9 m a Ilm - 1 = y/n m /\co\8 m <j"'"--\ witn c 
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' ) cr" (t) and rescaling the 



r 



with 



(2X//3h 2 — iXj/h) in the high temperature regime, leads to 
a more balanced EOM and can lead to a lower requirement in 
terms of the number of ADOs. The rescaled representation is 
useful for simulating large systems such as the FMO complex 
and for studying the ADOs themselves, as follows. 

To characterize the information flow from the point of view 
of the non-Markovian degrees of freedom, we propose another 
measure based on the trace distance. The auxiliary systems 
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FIG. 1. (a) Trace distance of the states p 1 = |1}(1| and p 2 = |2}(2| 
as a function of time for a two-site system (dimer) in the hierarchy 
equation of motion model. The trace distance evolves from maxi- 
mally distinguishable (= 1) for the initial states to indistinguishable 
(= 0) for the thermal equilibrium state. At intermediate times, the 
trace distance increases, which is due to a reversal of the informa- 
tion flow. Various bath correlation times are shown, 150 fs (blue), 
100 fs (red), and 50 fs (yellow). For long correlation times, that is, 
a bath with more memory, the regions of increasing trace distance 
are more pronounced, (b) The ADO distance Dado illustrates the 
information content in the non-Markovian degrees of freedom of this 
model for the same system, initial states, and bath correlation times 
as in (a). The anticorrelated oscillations of the two figures clearly 
show the information flow between the system and the NM degrees 
of freedom. 



describe perturbations on the time-evolved excitonic state due 
to the phonon environment. This separation into system and 
auxiliary degrees of freedom implies a quantification of the 
information content of the auxiliary systems at a given point 
in time: 



^ADO 



n^O 



(4) 



This measure depends on the electronic initial states and, im- 
portantly, uses the normalized ADOs. Although the <j n (t) 
are not quantum states for n / 0, the trace distance never- 
theless measures differences of these operators arising from 
the different initial states. This measure could in principle be 
extended to other families of master equations that employ 
auxiliary degrees of freedom. 

Results for a dimer— In this section, a two-molecule sys- 
tem, or dimer, is studied in terms of the measure for non- 
Markovianity, using, if not otherwise mentioned, the regu- 
lar hierarchy equation of motion approach Eq. (|3]l and scan- 
ning over the parameters of the model. As standard pa- 
rameters, we use the site energies e% — 0, £2 = 120/cm, 
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and the coupling J = — 87.7/cm. This corresponds to the 
strongly-coupled bacteriochlorophyll site 1 and 2 subsystem 
of the Fenna-Matthews-Olson complex. Along the lines of 
the discussion in lfl2l . the standard bath correlation time is 
taken from 50 fs to 150 fs. We assume ambient temperature 
T = 288 K (200/cm) . The main results use 40 tiers (819 
ADOs), and Eq. ([3]l is integrated up to maximally 20 ps. 

First, in Fig. [T] (a), we investigate the time evolution of 
the trace distance in Eq. ([T]i for the dimer system at reor- 
ganization energy A = 20/cm. We choose the two initial 
states p 1 = |1)(1| and p 2 = |2)(2| and the standard parame- 
ters given above, while varying the bath dissipation rate. One 
can see clear non-Markovian revivals in the time intervals at 
around 150 fs and around 300 fs, borne out by increases in 
the trace distance. For larger correlation time of the bath of 
t c = 150 fs these revivals are more pronounced. A bath ex- 
citation is more likely to return to the system instead of being 
dissipated away rapidly, as occurs in the case of smaller cor- 
relation time r c = 50 fs. 

In Fig. [T](b), we investigate the time evolution of the ADO 
distance measure in Eq. Q for the same system, using the 
normalized auxiliary systems. At small times, the ADOs are 
all zero, thus -Dado = 0> while at long times the system 
and ADOs converge to the thermal equilibrium and all infor- 
mation about the initial states is dissipated. At intermediate 
times, one finds strong information flow between system and 
the non-Markovian degrees of freedom, which is borne out by 
the anticorrelated oscillations in (a) and (b), respectively. For 
this, note the dips at around 150 fs and around 300 fs in (b). 
Additionally, the information content in the ADOs is smaller 
when the bath dissipation rate is larger. 

The computation of the non-Markovianity measure in Eq. 
Q involves an optimization of the initial states. On the one 
hand, the following results use the two initial states p x — 
|1)(1| and p 2 = |2)(2| as above, but one the other hand, 
compare to a systematically-optimized pair of initial states. 
One initial state is chosen out of 50 pure states on the Bloch 
sphere. This leads to a total number of 1225 independent pairs 
of states in the optimization. (The state optimization itself is 
performed with 20 tiers of ADOs.) 

We now proceed to scan the NM-ity measure over the pa- 
rameters of the model, beginning with the system parame- 
ters, the site energy difference and the chromophoric coupling. 
First, the NM-ity is strongly dependent on the chromophoric 
coupling, see Fig. [2] (a). At zero coupling, the molecules are 
independent and the NM-ity is zero; the two initial states re- 
main perfectly distinguishable. At non-zero coupling, the dy- 
namics shows increasing NM revivals, while both initial states 
converge to the same thermal equilibrium state. Second, with 
respect to the site energy difference, the dependence of the 
NM-ity weakly slopes downward, see Fig[2](b). As the en- 
ergy difference increases the coupling becomes less signifi- 
cant and the molecules become more independent. This leads 
to decreased NM-ity. 

Next, we investigate the role of the bath parameters. First, 
as mentioned before, the bath correlation time crucially af- 
fects the memory of the bath and the information back flow 
to the system. In Fig. [2] (c), we show the NM-ity mea- 



£ 0.35 

1 0.30 
> 0.25 
£ 0.20 
| 0.15 
f 0.10 
a 0.05 

2 0.00 





20 40 60 80 100120140 
Coupling [1/cm] 



50 100 150 200 250 300 
Site energy difference [1/cm] 





Reorganization energy [1/cm] 




20 40 60 80 100 120 140 
Bath dissipation rate [1/cm] 



20 40 60 80 100120 
Reorganization energy [1/cm] 



FIG. 2. Non-Markovianity as a function of the parameters of a dimer 
system, a)-d), and the Fenna-Matthews-Olson complex, e) and f), in 
the hierarchy equation of motion approach. For the dimer: a) cou- 
pling J12, b) site energy difference £2 — £i, c) bath dissipation rate 
7 (which is related to the bath correlation time by 7 = 1 /r c ,) and d) 
reorganization energy A. The solid lines are calculations performed 
with the initial states being p 1 — |1}(1| and p 2 — |2)(2| while the 
dots are calculations performed with optimized initial states. The 
dashed part of the lines indicates that a parameter was scanned be- 
yond the validity of the equation based on the truncation condition. 
For the FMO complex: e) bath dissipation rate and f) reorganization 
energy. In all figures, except in c) and e), we use the different bath 
correlation times: 50 fs (blue), 100 fs (purple), and 150 fs (red). 



sure as a function of the inverse bath correlation time, i.e. 
the bath dissipation rate. A longer r c leads to more non- 
Markovianity, because a sluggish bath keeps the memory of 
excitations much longer and the information is able to re- 
turn back to the electronic degrees of freedom. Second, the 
other significant bath parameter, the reorganization energy, is 
a measure of the non-equilibrium character of the phonon bath 
in the excited state. In Fig. |2](d), we show the NM measure 
as a function of the reorganization energy. At zero A, which 
essentially means no coupling to the phonons, the unitary dy- 
namics obviously shows zero non-Markovianity. At weak A, 
we observe a NM-ity of around 0.1 at t c = 150 fs, which is 
completely neglected by Redfield theory. At intermediate A of 
around 40/cm, close to the physiological values of the Fenna- 
Matthews-Olson complex, the NM-ity is maximal, showing 
a value of around 0.2. At large A, the regime of incoherent 
transport |12| . the NM-ity vanishes. 

Results for the Fenna-Matthews-Olson complex— We now 
investigate the exciton-phonon information flow in the Fenna- 
Matthews-Olson complex. The FMO complex is found in 
green sulphur bacteria, where it connects the antenna complex 
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with the reaction center during the energy transfer process. 
The FMO complex is a trimer with seven bacteriochlorophyll 
(BChl) molecules in each subunit and with one additional 
BChl molecule shared between the units. We implement the 
hierarchy equation of motion for the seven sites (N = 7), 
using the method for rescaling the auxilliary systems as dis- 
cussed above and in [16|. As was shown recently ifTTl . with 
four tiers (329 ADOs) the coherent population dynamics can 
be accurately simulated, but due to the reduced number of 
ADOs, we expect to minimally overestimate the NM-ity. We 
use the initial states p 1 = |1)(1| and p 2 = |2)(2|. First, the 
scan of 7 obtains a strong dependence similar to the dimer 
system, see Fig. |2](e). Second, the scan of the reorganization 
energy obtains a maximum at around A = 55/cm with the 
NM-ity being 0.2 at A = 35/cm and r c = 150 fs, see Fig. 
|2](f). The results suggest that non-Markovian effects should 
play a significant role in the function of the Fenna-Matthews- 
Olson complex. 

Conclusion— The theoretical and experimental characteri- 
zation of photosynthetic light-harvesting complexes in terms 
of non-Markovianity, as e.g. defined in the recent work 



lfT3l[T4ll . can lead to a greater understanding of the exciton dy- 
namics in these systems. In principle, the recently proposed 
ultrafast quantum process tomography protocol extracts a full 
characterization of the quantum map lfl"8l . and could thus be 
used for the experimental side of this task. However, experi- 
ments specifically designed to extract a particular observable, 
such as in our case the NM-ity measure, could be carried out 
with a reduced number of experimental requirements. Future 
work will investigate analytically and numerically how to effi- 
ciently extract the non-Markovianity from the phase-matched 
signal observed in four-wave mixing experiments. 

To summarize, we have shown numerically that there is a 
considerable non-Markovian information flow between elec- 
tronic and phononic degrees of freedom of light-harvesting 
chromophoric systems under physiological conditions, utiliz- 
ing a realistic model for systems such as the Fenna-Matthews- 
Olson complex. The results suggest that non-Markovianity 
plays a significant role in photosynthetic energy transfer. 
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